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ABSTRACT 

In  this  paper  we  study  strictly  positive  definite  functions  on  the  unit  sphere  of  the 
m-dimensional  Euclidean  space.  Such  functions  can  be  used  for  solving  a  scattered  data 
interpolation  problem  on  spheres.  Since  positive  definite  functions  on  the  sphere  were 
already  characterized  by  Schoenberg  some  fifty  years  ago,  the  issue  here  is  to  determine 
wlmt  kind  of  positive  definite  functions  arc  actually  strictly  positive  definite.  The  study 
of  this  problem  was  initiated  recently  by  Xu  and  Cheney,  [XC],  where  certain  sufficient 
conditions  were  derived.  A  new  approach,  which  is  based  on  a  critical  connection  between 
this  problem  and  that  of  multivariate  polynomial  interpolation  on  spheres,  is  presented 
here.  The  relevant  interpolation  problem  is  subsequently  analysed  by  three  different  com¬ 
plementary  methods.  The  first  is  based  on  the  de  Boor-Ron  general  “least  solution  for 
the  multivariate  polynomial  interpolation  problem”.  The  second,  which  is  suitable  only 
for  m  =  2,  is  based  on  the  connection  between  bivariate  harmonic  polynomials  and  uni¬ 
variate  analytic  polynomials,  and  reduces  the  problem  to  the  structure  of  the  integer  zeros 
of  bounded  univariate  exponentials.  Finally,  the  last  method  invokes  the  realization  of 
harmonic  polynomials  as  the  polynomial  kernel  of  the  laplacian,  thereby  exploiting  some 
basic  relations  between  homogeneous  ideals  and  their  polynomial  kernels. 
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Strictly  Positive  Definite  Functions  on  Spheres 

Amos  Ron  &  Xingping  Sun 

1.  Introduction 

Let  5m_1  denote  the  unit  sphere  in  the  Euclidean  space  IRm  (m  >  2).  and  drn  the 
geodesic  distance  on  Sm~l,  i.e., 

dm(x,y)  =  Arccos(a;  •  y),  x.j/tS"1'1. 

Here  x  ■  y  denotes  the  usual  inner  product  of  x  and  y.  Let  g  :  [0, 7r]  — >  ]R  be  a  continuous 
function,  and  let  0  =  {0lr . . . ,  9n}  be  a  subset  of  Sm~l.  We  study  in  this  paper  the  possible 
strict  positive  definiteness  of  the  n  x  n  matrix 

(1.1)  A  :=  A©  :=  Ag<& 
with  ij —entry 

(1.2)  Aij  =  g(dm(0uej)); 

i.e.,  we  look  for  conditions  such  that 

n  n 

cT Ac  =  Y^'Y^cici9(dm{Qi,Qj))  >0,  Vc  =  (ci,...,cn)  6  !Rn\0. 

i=i  j=i 

The  matrix  A  of  (1.1)  naturally  arises  in  the  study  of  approximations  to  scattered 
data  on  spheres.  Let  ©  C  5m_1  be  a  finite  set,  and  let  /  be  some  function  defined  on  0. 
To  approximate  /  by  a  function  defined  on  the  entire  sphere,  one  may  choose  a  univariate 
function  g  :  [0, 7r],  and  look  for  an  interpolant  gj  in  the  linear  space 

G©  :=  span{s(dm(-,0))  :  0  6  0}. 

The  existence  of  a  unique  interpolant  gj  6  G©  for  /  then  amounts  to  the  invertibility  of 
the  above  matrix  A.  Of  course,  if  A  is  also  positive  definite,  then  the  finding  of  gj.  i.e., 
the  inversion  of  A,  can  then  be  approached  by  efficient  and  stable  numerical  (iterative  or 
direct)  methods. 

In  view  of  the  above,  the  following  problem  becomes  self-suggestive: 


Problem  1.3.  Determine  conditions  under  which  the  interpolation  matrix  .4©  of  (1.1)  is 

(a)  Positive  definite  (for  Sm_1),  for  any  ©  C  Sm-1  of  cardinality  n,  for  some  fixed  n. 

(b)  Positive  definite  (for  for  any  ©  C  Sm~1. 

(c,d)  Same  as  (a ,b),  with  “ invertibility ”  replacing  “ positive  definiteness”. 

Definition  1.4.  Let  g  be  a  univariate  continuous  function  defined  on  [0,  ir\.  U'e  say  that  g 
is  (strictly)  positive  definite  of  order  n  for  Sm~l  if  for  each  ©  C  5m_1  of  cardinality 
n  the  corresponding  matrix  .4©  is  ( strictly )  positive  definite.  A  function  g,  that  is  (strictly) 
positive  definite  of  all  orders,  is  (strictly)  positive  definite. 

The  analogous  problem  in  Euclidean  spaces,  i.e.,  when  ©  C  IRm,  has  been  intensively 
studied  in  the  literature.  In  [Si],  Schoenberg  had  characterized  the  positive  definite  func¬ 
tions  of  all  orders  for  IRm,  and  in  [M].  Micchelli  established  the  invertibility  of  certain 
interpolation  matrices  arising  from  approximating  scattered  data  in  IRm.  Micchellis  re¬ 
sults  had  led  to  a  wealth  of  results  (cf.  e.g.  [NW1,2],  [QSW],  [SiWa]),  in  which  estimates 
for  various  norms  and  corresponding  condition  numbers  of  the  interpolation  matrix  .4  were 
established.  Wo  refer  to  the  review  article  of  Dyn  [D]  for  more  information  on  this  subject. 
Schoenberg  also  considered  the  problem  on  the  sphere.  In  [S2]  he  proved  the  following 
characterization  for  a  function  g  to  be  positive  definite  of  all  orders.  In  this  result.  P ^ 
denotes  the  kth  degree  Gegenbaucr  (“ultraspherical”)  polynomial  associated  with  A.  [Sz, 
p.81],  [SW,  p.  148]. 

Result  1.5  [S2].  A  continuous  function  g  :  [0, 7r]  is  positive  definite  on  Sm~l  if  and  only 
if  it  has  the  form 

00 

(1.6)  g(t)  =  ^2akP(kx){cost), 

k= 0 


in  which  A  =  (m  —  2)/2,  a*  >  0,  and  akP^\  1)  <  oo. 

Our  interest  in  the  problem  was  initiated  by  the  paper  [XC]  of  Xu  and  Cheney,  in 
which  the  following  question  is  addressed:  “find  conditions  on  the  coefficients  (afc)feeN  in 
(1.6),  under  which  g  is  strictly  positive  definite”  (either  of  a  specific  order  or  of  all  orders). 
Among  various  other  results,  it  is  shown  in  [XC]  that,  if  all  the  coefficients  ak  in  (1.6)  are 
positive,  then  the  function  g  is  strictly  positive  definite  on  Sm_1.  Further  discussions  of 
this  problem  can  be  found  in  [LC]  and  [CX]. 

As  this  paper  will  show,  a  close  relationship  exists  between  the  problem  of  determining 
strict  positive  definiteness  and  that  of  multivariate  polynomial  interpolation.  This  connec¬ 
tion  is  discussed  in  §2,  and  allows  us  to  find  an  equivalent  version  to  the  original  problem 
in  terms  of  polynomial  interpolation  on  the  unit  sphere. 
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We  then  present  three  different  methods  for  analysing  the  equivalent  polynomial  in¬ 
terpolation  problem.  The  first  exploits  the  de  Boor-Ron  “least  solution  for  the  polynomial 
interpolation  problem"  (cf.  [BR1-3]).  This  method  yields  that  g  is  strictly  positive  definite 
of  order  n  if  (but  not  only  if)  the  corresponding  coefficients  ao,a\. a.[n/2\  are  positive 
(cf.  Theorem  3.2.  The  proof,  as  well  as  a  complementary  discussion,  are  provided  in  §4). 
For  the  case  m  =  2,  this  result  was  already  proved  in  [XC].  Another  approach  takes  advan¬ 
tage  of  the  connection  between  spherical  harmonics  on  the  circle  and  analytic  polynomials, 
and  allows  us  to  characterize,  for  m  =  2.  the  strict  positive  definiteness  of  g  in  terms  of  the 
positive  integer  zeros  of  univariate  exponentials  with  real  coefficients.  Our  principal  find¬ 
ings  that  results  from  this  approach  are  collected  in  Theorem  3.5.  with  proofs  and  further 
discussion  to  be  found  in  §5.  Finally,  in  §6.  we  choose  a  different  track  which  is  suitable  for 
any  spatial  dimension,  and  makes  use  of  the  fact  that  the  harmonic  polynomials  are  the 
polynomials  in  the  kernel  of  the  laplacian.  That  direction  utilizes  some  of  the  polynomial 
ideal  basics,  and  its  main  results  are  collected  in  Theorems  3.6  and  3.7.  These  theorems 
imply,  in  particular,  that  g  is  strictly  positive  definite  if  the  set  I\  :=  {k  6  7Z+  :  a k  >  0} 
(with  (ak)  as  in  1.6)  contains  arbitrarily  long  sequences  of  consecutive  even  integers,  as 
well  as  arbitrarily  long  sequences  of  consecutive  odd  integers. 

We  use  in  the  paper  the  following  notations  and  conventions. 

n 

stands  for  the  space  of  all  polynomials  in  m  variables  (where  the  value  of  m  should  be 
clear  from  the  context).  The  subspace  of  II  that  consists  of  the  fc-degree  homogeneous 
polynomials  is  denoted  by 

nj. 

Also,  given  any  I(  C  ZZ+,  we  set 

m-  :=  £  nj- 

k£K 

If  K  =  {0, 1, 2, ... ,  k}.  we  often  use  11*  :=  lift-,  i.e.,  life  is  the  space  of  all  polynomials  of 
degree  <  k.  A  parallel  set  of  notations  is  used  for  harmonic  polynomials.  Here,  we  set 

n 

for  the  space  of  all  harmonic  polynomials,  and  define 

n°k  -.=  n£  n  n,  nK  ■=  n K  n  n,  nk  ■.=  nk  n  n. 

Filially,  an  exponential  in  this  paper  is  either  a  function  of  the  form 
(1.7)  eg  :JRm  ->C  e9x .  0  €  Cm, 
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or  any  finite  linear  combination  of  such  functions. 

Additional  definitions  are  used  in  the  material  presented  in  §6,  and,  to  a  lesser  extent, 
in  the  last  part  of  §4.  Most  of  these  are  given  in  §3  (prior  to  the  statement  of  Theorem 
3.6,  which  is  taken  from  §6). 


2.  The  connection  to  polynomial  interpolation  by  harmonic  polynomials 

Since  we  are  seeking  conditions  that  characterize  the  strict  positive  definiteness  of 
g ,  and  since  Schoenberg’s  theorem  already  characterizes  the  positive  definiteness  of  g.  we 
may  assume  that  g  has  the  form  (1.6).  In  turn,  that  allows  us  to  write  the  matrix  .4©  as 
the  infinite  sum 


(2.1) 


oo 

A©  =  ^2  ak^k, 

k=0 


with 

A,.:=  (PW(0.*))M€0. 

Since  each  summand  A k  is  positive  definite  (by  virtue  of  Schoenberg’s  result,  but  also  in 
view  of  (2.3)  below),  this  immediately  shows  that  the  matrix  A  is  more  likely  to  be  strictly 
positive  definite,  with  the  increase  of  the  non-zero  coefficients  (at)fc  in  the  representation 
(1.6)  of  g.  Also,  since  the  matrix  A k  is  positive  definite,  it  can  be  written  in  the  form 
CjCk-  Among  such  factorizations  of  A) t,  we  select  below  a  particular  one,  which  is  based 
basic  properties  of  spherical  harmonics,  and  which  will  allow  us  to  draw  the  connection 
between  Problem  1.3  and  polynomial  interpolation  on  the  sphere. 

Recall  that  the  restriction  to  5m_1  of  a  homogeneous  harmonic  polynomial  of  degree  k 
is  called  a  spherical  harmonic  of  degree  k.  The  following  fact  about  spherical  harmonics 
can  be  found,  for  example,  in  Stein  and  Weiss  [SW,  Chapter  IV]. 

Let  . . . ,  }  be  an  orthonormal  basis  ofH°k  (here,  hk  :=  dimH,®).  Then,  there 

is  a  positive  constant  Ck,\  such  that 


(2.2)  P?\x  ■  y)  =  >(!,). 

j  =  1 

Here  X  =  (rn  —  2)/2,  and,  as  before,  P ^  is  the  appropriate  Gegenbauer  polynomial. 
Using  (2.2),  we  observe  that  Ak  of  (2.1)  can  be  factored  as  follows: 

(2.3)  Ak  =  Ck,\BlBk, 
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with 

(2.4)  Bk  :=  Bk,e  :=  (c,e. 

Therefore,  for  any  vector  c  e  1R  ,  cTAkC  =  ||Bfcc||^,  and  hence 

oo 

(2.5)  ctAc  =  ak\\Bkc\\l. 

k= 0 

This  proves  one  of  the  implications  in  Schoenberg’s  theorem,  and  implies  that  the  problem 
of  whether  g  is  strictly  positive  definite  or  not  depends  only  on  the  set 

(2.6)  Km,g  :=  {k  e7Z+  :  ak  ±  0}. 
and  not  on  the  actual  (positive)  value  of  each  ak,  k  £  Km%g. 

Proposition  2.7.  Let  g\  and  g>  be  two  positive  definite  functions  on  If  Km.gi  = 

I\,n,g2  •  then,  for  every  0  C  STn~l,  the  matrix  Agiie  is  invertible  if  and  only  if  the  matrix 
Ag  ,  ,©  is  invertible. 

Because  of  the  above  proposition,  we  modify  our  Problem  1.3  as  follows: 

Modified  Problem  2.8.  Study  the  following  notion:  Given  a  subset  K  C  Z+.  we  say 
that  I\  induces  strict  positive  definiteness  (of  order  n)  for  Sm~1  if  every /some 
positive  definite  g  that  satisfies  Km,g  =  K  is  strictly  positive  definite  (of  order  n). 

Remark  2.9.  Note  that  whenever  g  is  strictly  positive  definite  for  Sm~l,  it  is  certainly 
so  for  Sl~l,  l  <  m.  In  view  of  Schoenberg’s  result,  it  can  therefore  be  expressed  as  a 
linear  combination  of  Gegenbauer  polynomials  of  smaller  orders.  In  this  regard,  it  is  useful 
to  recall  the  following  connection  between  Gegenbauer  polynomials  of  different  orders  (cf. 
[AF]):  for  every  A  >  v  >  0,  and  a  non- negative  k,  there  exist  positive  coefficients  (a\Mikjn)n 
such  that 

pW  -  v  *  p{,/) 

0<2 n<k 

As  stated  before,  we  begin  our  study  of  Problem  2.8  by  linking  it  to  the  problem  of 
interpolating  by  spherical  harmonics  on  the  sphere. 

Theorem  2.10.  Given  K  C  Z Z+,  and  a  positive  definite  g  that  satisfies  Km<g  -  K.  the 
following  conditions  are  equivalent,  for  any  0  C  Sm~l. 

(a)  The  matrix  A9iq  is  invertible. 
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(b)  The  restriction  oflin  to  ©  is  of  full  dimension  #©,  i.e.,  every  f  :  0  can  be  interpolated 
by  a  polynomial  p  6  Hk- 

Proof.  By  (2.5),  A9<q  is  singular  if  and  only  if  BkC  =  0,  for  some  c  €  IRe\0,  and 
for  all  k  €  K.  This  is  equivalent  to  the  saying  that  the  linear  functional  <p  :=  Qco$e 
(with  the  evaluation  at  6)  vanishes  on  a  basis  for  each  TL°k,  k  e  K,  hence  vanishes  on 
Hk-  Finally,  the  existence  of  a  non-zero  linear  functional  that  is  supported  on  0  and 

annihilates  Hk  is  equivalent  to  the  inequality  dim('H^'|e)  <  #0.  4 

3.  Main  Results 

Lemma  3.1  is  proved  in  §4.  Theorem  3.5  follows  from  Theorem  5.1,  Corollary  5.2, 
Theorem  5.3,  and  Corollary  5.7  of  §5.  Theorems  3.6  and  3.7  are  proved  in  §6.  No  formal 
proof  is  required  for  the  other  results  stated  in  this  section. 

The  space  n©  mentioned  in  following  statement  is  the  least  solution  of  the  polynomial 
interpolation  problem,  [BRl]. 

Lemma  3.1.  For  any  0  C  5”*1  of  cardinality  n,  the  following  is  true: 

(a)  There  exists  a  space  n©  C  H[n/2}  such  that  interpolation  from  II©  to  any  f  defined 
on  0  is  always  possible,  i.e.,  the  restriction  of  II©  to  0  is  of  dimension  n  =  #0. 

(b)  If  n  is  even,  and  0  does  not  lie  on  a  circle,  then  II©  C  Hn/2-i  ~  %|n/2J-i- 

(c)  If  n  ^  5  is  odd  and  no  (n  —  1)  points  of  ©  lie  on  a  circle,  then  II©  C  7f(n-3)/2  = 

^  Ln/2J  - 1  ■ 

Combining  Theorem  2.10  with  (a)  of  Lemma  3.1,  we  arrive  at  the  following  result. 

Theorem  3.2.  Let  K  C  7L+.  Then,  given  any  m  >  2,  K  induces  strict  positive  definite¬ 
ness  of  order  n  for  Sm~1  if  {0, 1, ... ,  \  n/ 2J }  C  K. 

For  m  —  2,  the  above  theorem  is  due  to  [XC].  Note  that  the  theorem  reproduces  the 
Xu-Cheney  result  which  asserts  that  g  is  strictly  definite  positive  whenever  Kmtg  =  2Z+. 
The  statement  of  the  theorem  is  sharp  in  the  sense  that  the  set  K  :=  {0, 1, . . . ,  [n/2j  —  1} 
does  not  induce  strict  positive  definiteness  of  order  n.  However,  parts  (b)  and  (c)  in 
Lemma  3.1  imply  that,  for  m  >  2,  this  latter  K  “almost”  induces  the  required  strict 
positive  definiteness. 

Corollary  3.3.  Let  g  be  a  positive  definite  function  for  5m_1,  K  :=  Km,g.  Assume  that 
{0. 1, . . . ,  [n/2j  -  1}  C  K,  and  let  ©  C  5m_1  be  of  cardinality  n.  Then  the  corresponding 
matrix  A©  of  (1.1)  is  strictly  positive  definite  if  one  of  the  following  conditions  holds: 

(a)  n  is  even,  ©  does  not  lie  on  a  circle. 


(b)  n  is  odd,  and  no  n  -  1  points  from  ©  lie  on  a  circle. 

Applying  Remark  2.9  to  Theorem  3.2,  we  arrive  at  the  following  corollary.  A  different 
proof  for  that  corollary,  that  follows  from  the  results  of  §6,  is  sketched  at  the  *>nd  of  the 
present  section. 

Corollary  3.4.  If  g  is  positive  definite  for  Sm~l,  and  if  Km<g  contains  one  even  integer 
k i,  as  well  as  one  odd  integer  kj,  then  g  is  strictly  positive  definite  of  order  n  for  Sl~l, 
provided  that  l  <  m,  and  [n/2j  <  min(Aq,  A^). 

Theorem  3.2  is  useful  if  one  wants  to  assert  the  strict  positive  definiteness  of  g  of  a 
certain  order  n,  since  it  allows  g  to  be  a  polynomial  of  small  degree.  However,  the  theorem 
provides  no  real  clue  to  the  problem  of  strict  positive  definiteness  (of  all  orders)  of  g ,  since 
it  only  implies  (the  aforementioned  result)  that  g  is  strictly  positive  definite  whenever 
I\m,g  =  TL+.  Results  concerning  the  induction  of  strict  positive  definiteness  of  all  orders 
are  derived  in  §5  and  §6. 

Section  5  is  devoted  to  the  case  m  —  2.  Some  of  its  findings  are  collected  in  our  next 
theorem. 

Theorem  3.5.  Let  K  C  ZZ+.  Then, 

(a)  I\  induces  strict  positive  definiteness  for  Sl  if  and  only  if  there  exists  no  non-zero 
function  of  the  form 

£>eir,  T  C  (0,27r),  #T  <  oo,  cT  €  IR.  Vr 
re  T 

that  vanishes  on  K. 

(b)  In  order  that  I\  C  ZZ+  induce  strict  positive  definiteness  for  Sl ,  it  is  necessary  that 
K  has  infinite  intersection  with  any  set  of  the  form  k7Z+,  k  G  IX.  The  same  applies  to 
sets  of  the  form  k/2  +  k7L+  (provided  that  k  is  even).  However,  I\  need  not  intersect 
at  all  sets  of  the  form  a-\-k7L+,  if  a  ^  0,  kj  2.  In  fact,  for  such  a,  I\  :=  7Z+\(a  +  k7Z+ ) 
induces  strict  positive  definiteness  for  Sl. 

(c)  In  order  that  K  C  ZZ+  induce  strict  positive  for  S1,  it  suffices  that,  for  every  finite 
set  J  in  IN\{1},  and  any  n  G  IN  there  exist  a  G  ZZ+  and  k  G  1N\(J1X),  such  that 
{a,  a  4-  k, . . . ,  a  +  (n  —  l)fc}  C  K. 

In  §6,  we  investigate  the  case  of  general  m.  There,  we  use  the  following  additional 
notations  and  terminology. 

The  first  required  notion  is  that  of  exponential  spaces.  These  are  defined  in  terms 
of  some  set  Cl  C  lRm,  and  a  positive  integer  n  as  follows: 

En(f2)  :=  {£  c0ee  :  ©  C  ft,  #©  <  n,  ce  €  C,  V0}. 
eee 
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(Note  that  eg  is  a  constant,  but  eg  is  a  function.)  Obviously,  the  above  exponential  spaces 
are  not  linear  spaces,  in  contrast  with  the  larger  space 


E{Q)  :=Un>i En(Sl). 

Another  set  of  notations  concerns  maps  defined  on  the  algebra  A  of  all  formal  power 
series  in  m-variables.  Since  A  is  the  direct  (infinite)  sum  J2k> o  ^2  °f  the  homogeneous 
polynomial  spaces,  there  exists,  for  every  k  >  0,  a  well-defined  projector 

k:  A  -»  ng, 

that  assigns  to  each  power  series  /  €  A  the  fc-degree  homogeneous  component  in  its  power 
expansion.  Further,  given  a  subset  K  C  Z+,  the  sum 

K:=^2~k 

k£l< 

defines  an  analogous  projector,  this  time  from  A  onto  FI /v- . 

Finally,  we  reserve  a  special  notation  T  for  the  polynomial 

m 

rir>  :=  rm(i)  :=  v  4. 

n=  1 

Thus,  TA  is  the  (homogeneous  principal)  ideal  generated  by  T.  Note  that  the  evaluation 
r (D)  of  r  at  D  is  the  laplacian  operator. 

The  principal  result  of  §6  is  as  follows: 

Theorem  3.6.  Let  I\  C  5Z+,  and  let  n  be  a  positive  integer.  Then  K  does  not  induce 
strict  positive  definiteness  of  order  n  for  Sm~l,  if  and  only  if  there  exists  f  E  En(Sm~l) 
for  which  K(f)  G  TA. 

In  fact,  the  proof  of  Theorem  3.6  is  more  informative  than  its  statement:  given  ©  C 
STn~i,  the  proof  shows  that  the  matrix  A9t@  is  singular  if  and  only  if  there  exists  an 
exponential  /  G  E(Q)  that  satisfies  Km,g(f)  G  T A. 

The  following  sufficient  condition  is  derived  in  §6  from  the  characterization  in  Theorem 
3.6. 

Theorem  3.7.  Given  0  C  Sm~l,  define 

cr(0)  :=  min{#©'  C  ©  :  span(0\0')  £  IRm}. 

Let  j  be  the  minimal  integer  that  satisfies  >  <x(©)-  Let  g  be  a  positive  defi¬ 
nite  function  for  Sm~l ,  and  assume  that  the  set  { k  G  :  k  >  (#©)/ 2}  contains  j 

consecutive  even  integers  as  well  as  j  consecutive  odd  integers.  Then  .49,©  is  invertible. 

Since  cr(0)  <  n  —  m  +  1  for  any  ©  C  lRm  of  cardinality  n,  we  obtain  the  following 
corollary: 
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Corollary  3.8.  Let  K  C  7Z+.  and  n  a  positive  integer.  Then  K  induces  strict  positive 
definiteness  of  order  n  on  Sm~l,  if,  with  j  the  minimal  integer  that  satisfies  (  jr_;2)  > 
n-m  +  1,  there  are  j  consecutive  even  integers  and  j  consecutive  odd  integers  in  the  set 
{A:  >  nj 2}  H  K. 

In  particular,  if  K  contains  arbitrarily  long  sequences  of  consecutive  evens  and  of 
consecutive  odds,  then  K  induces  strict  positive  definiteness  (of  all  orders)  on  Sm~1,  for 
every  m  >  2. 

Theorem  3.7  also  provides  another  proof  for  Corollary  3.4:  if  ©  lies  on  Sl~\  l  <  m. 
then,  in  particular,  it  lies  on  linear  proper  subspace  of  IRm.  This  implies  that  cr(Q)  =  0, 
and  one  can  take  in  Theorem  3.7  j  =  1. 


4.  The  least  solution  of  the  polynomial  interpolation  problem 

We  prove  here  Lemma  3.1. 

Given  any  polynomial  of  degree  k,  we  denote  by  the  unique  polynomial  in  11°.  that 
satisfies  deg(p  —  pj)  <  lc,  and  refer  to  pj  as  the  leading  term  of  p.  Also,  given  p  €  II.  we 
use  p{D)  to  denote  the  corresponding  constant  coefficient  differential  operator.  Directional 
derivatives  are  denoted  by  Dx.  where  x  €  IRm  is  the  direction. 

Given  any  finite  0  C  HTn,  the  paper  [BR1]  introduces  a  polynomial  space  n©  that 
satisfies  all. the  following  properties: 

(a)  For  any  function  /  :  0,  there  exists  a  unique  p  €  II©  that  agrees  with  /  on  0. 

(b)  II©  is  D-invariant.  i.e.,  for  every  x  €  Dx II©  C  II©  (equivalently,  II©  is 

translation-invariant ) . 

(c)  II©  is  homogeneous.  That  is 

OO 

n©  =  ne>J, 

j=0 

where  11©^  :=  II©  H  11°. 

(d)  If  the  polynomial  p  vanishes  on  0,  then  its  leading  term  pj  annihilates  II©  in  the  sense 
that  pf  (Z?)II©  =  0. 

(e)  Conversely,  every  homogeneous  polynomial  that  annihilates  II©  is  the  leading  term  of 
some  other  polynomial  that  vanishes  on  0. 

We  refer  to  [BRl]  and  [BR3]  for  more  details  about  II©.  The  exact  definition  of  II© 
will  be  given  in  the  sequel,  and  will  be  used  for  the  proof  of  (c)  in  Lemma  3.1.  For  the 
time  being,  though,  we  need  only  the  fact  that  a  polynomial  space  that  satisfies  these  five 
properties  exists. 
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Lemma  4.1.  Let  k  he  the  least  integer  that  satisfies  lie  C  11*.  If  0  C  S'11-1,  then 

dim  II©  j  >  2.  j  -  1,2,  1. 

Claim  (a)  of  Lemma  3.1  follows  from  Lemma  4.1.  Indeed,  since  n©  is  translation- 
invariant,  it  must  contain  the  constants,  hence  dimll©^  =  1.  Also,  by  the  definition  of  k, 
dimn©tA:  >  1.  Hence,  assuming  Lemma  4.1  has  been  proved,  we  obtain  that 

#0  =  dimll©  >  1  +  2(Jb  -  1)  +  1  =  2k, 

with  the  right-most  equality  due  to  property  (a)  of  II©. 

Proof  of  Lemma  2.2.  We  assume  that  the  claim  of  the  lemma  is  false,  and  will  seek 
a  contradiction.  Let  0  <  j  <  k  be  the  maximal  integer  that  violates  the  lemma's  claim.  If 
j  =  k  —  1,  then  IIqj+i  -  n©.*,  ^  0,  by  the  definition  of  k.  Otherwise,  by  the  maximally 
of  j,  dimn©,j  +  !  >  2.  Either  way.  n©,J  +  i  contains  a  non-zero  polynomial  p.  We  consider 
the  map  tp  :  IRm  — >  n”  defined  by 

V’(.r)  =  Dxp. 

Since  LI©  is  D- invariant,  rani'  C  II©,  hence  ran 0  C  n©j.  Since  degp  >  0,  0  ^  0.  On  the 
other  hand,  dim  ran  i/;  <  dimn©.j  <  1.  Therefore,  dim  ran0  =  1,  hence  dimker  0  =  m  -  1. 
Since  0  is  linear,  its  kernel,  then,  is  a  hyperplane,  and  p  has  to  be  a  univariate  polynomial 
of  the  form 

p(y)  =  c(Z-yy+1, 


with  £  perpendicular  to  the  above  hyperplane. 

Let  T (D)  be  the  m-dimensional  Laplacian.  We  note  that,  since  j  >  1,  T{D)p  = 
d  JliLi  tiit '  y)j~l  ^  On  t^ie  other  hand,  the  quadratic  polynomial  q  :=  T  -  1  vanishes 
on  0,  hence  by  property  (d)  of  II©,  T(D)  ~  q\{D)  annihilates  the  entire  n©.  In  particular. 
T(D)p  =  0,  which  contradicts  the  previous  conclusion.  4* 

The  proof  of  (b)  in  Lemma  3.1  is  also  quite  simple.  We  present,  however,  a  slightly 
longer  proof,  which  prepares  also  for  the  proof  of  (c).  We  assume  that  n  is  even,  and 
that  II©  £  nri/2_i,  and  will  prove  that  0  lies  on  a  circle.  First,  by  (a)  of  Lemma  3.1, 
II©  C  nn/2.  Second,  we  have  dimll©^  =  1,  and,  since  we  assume  that  dim!!©  „/2  >  1, 
Lemma  4.1  implies  that  dimllej  >  2,  j  =  1, 2, . .  .,n/2  -  1.  Taking  into  account  the  fact 
that  n©  is  n-dimensional,  we  realize  that  none  of  these  inequalities  can  be  sharp,  i.e.,  the 
homogeneous  dimensions  of  FI©  are 


1,2,2,. ..,2,1. 
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Let  q  €  IIe,j\0,  1  <  j  <  n/2,  and  let  tpq  be  the  corresponding  map  that  was  introduced 
in  the  proof  of  Lemma  4.1.  The  proof  of  that  lemma  shows  that  ran ipq  c  Usj-i-  Since 
dimllej-x  =  2,  this  implies  that  dimker ipq  >  m- 2,  i.e.,  th=»t  q  is  a  bivariate  polynomial. 
Further,  the  argument  in  Lemma  4.1  makes  clear  that  q  cannot  be  univariate.  Conse¬ 
quently,  rank  tpq  =  2,  and  ran  wq  =  Ilej-i-  By  selecting  p  to  be  any  non-zero  polynomial 
in  n0  n/2,  and  choosing  a  €  ZZ+  such  that  |a|  =  n/2  —  j  and  q  :=  Dap  ^  0,  we  obtain  that 
each  n©,_i,  j  =  2, . . .  ,n/ 2.  is  generated  by  derivatives  of  p.  That  same  result  is  trivial 
for  j  =  n/ 2  +  1  (since  n©iU/2  is  1-dimensional),  and  for  j  =  1.  Thus, 

n©  =  {P(D)P  :Pe  n}, 

i.e..  n©  comprises  of  the  derivatives  of  the  single  polynomial  p,  hence,  in  particular,  II© 
is  annihilated  by  the  (m  —  2)-dimensional  space  {Dx  :  x  £  ker^'p}.  Selecting  any  basis 
for  that  space,  we  invoke  property  (e)  of  II©  to  conclude  that  there  exist  m  —  2  linearly 
independent  linear  polynomials,  each  of  which  vanishes  on  ©,  i.e..  that  0  lies  on  a  2- 
dimensional  linear  manifold.  Since  ©  is  assumed  to  lie  also  on  STn~ 1 ,  we  conclude  that, 
indeed,  it  lies  on  a  circle.  This  completes  the  proof  of  (b). 

Before  we  prove  (c),  we  mention  that  its  statement  is  sharp  in  the  following  sense. 
First,  if  0  consists  of  five  points  in  S2,  then,  regardless  of  its  distribution,  the  5-dimensional 
n©  cannot  lie  in  the  4-dimensional  II(n_3)/2  =  III,  i.e.,  the  statement  in  (c)  fails  to  hold 
for  n  =  5.  Further,  if  all  the  points  of  0  except  one  lie  on  a  circle,  then  II©  qL  n(n_3)/2, 
and,  therefore,  one  cannot  prove  in  (c)  that  ©  must  entirely  lie  on  a  circle. 

In  order  to  prove  (c),  we  assume  that  n  >  7  is  odd,  and  that  II©  <£  n(n_3)/2  (the 
claim  is  trivial  for  n  =  1,3).  We  need  to  prove  then  that,  save  perhaps  one  point,  0  lies 
on  a  circle.  Here,  we  need  to  recall  from  [BR1]  that  the  definition  of  n©  is 

(4.2)  II©  =  spanf/i  :  /  e  £(©)}, 

with  E(Q)  (as  before)  the  span  of  the  exponentials  9  6  0,  and  with  f[  the  smallest 
degree  non-zero  homogeneous  term  is  the  power  expansion  of  /,  (i.e.,  with  k  the  minimal 
non-negative  intege  •  ^hat  satisfies  k(f)  #  0,  fl  :=  &(/)).  It  is  then  easy  to  check  that,  if  0 
contains  four  points  which  are  not  co-planar,  then  dimll©,!  >  3.  Therefore,  assuming  0 
not  to  lie  on  a  circle,  we  must  have  that  dimn©,i  >  3.  Repeating  the  same  counting  argu¬ 
ments  that  we  employed  in  the  proof  of  (b),  we  conclude  that  the  homogeneous  dimensions 
of  II©  must  be 

1,3, 2, 2,. ..,2,1. 

Selecting  any  p  €  II© ,(n_!)/2,  and  repeating  the  argument  that  was  used  in  the  proof  of 
(b),  we  conclude  that  p  is  bivariate  and  that  the  derivatives  of  p  form  a  subspace  in  n©  of 
dimension  n  —  1  (the  argument  relies  on  the  fact  that  dimn©  (n_1)/2-i  =  2,  hence  requires 
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n  >  7).  Let  x  €  IRm  be  such  that  Dxp  =  0.  Let  also  /  €  £(©)  be  an  exponential  that 
satisfies  fl  =  p.  Set  g  :=  Dxf ,  k  :=  (n  -  l)/2.  Since  j(/)  =  0,  for  j  <  k.  we  have  that 
U  ~  1  )(g)  =  0.  Also,  since  k(f)  =  p,  we  have  (k  -  1  )(g)  =  (k  -  1  ){Dxf)  =  Dxk{f)  = 
ArP  =  0.  Therefore,  deg(pt)  >  A:.  Since  g j  is  in  lie,  and  lie  C  II*,  we  must  thus  have 
deg(pi)  =  k.  hence  that  is  a  constant  multiple  of  p  (since  dim  II©,*  =  1),  say,  g t  =  cp. 
This  implies  that,  with  q  the  linear  polynomial  q(y)  :=  x  ■  y  -  c,  j{q(D)f )  =  0,  0  <  j  <  k, 
and  thereby  that,  if  q(D)f  ±  0,  then  deg (q(D)f)[  >  k,  in  contradiction  with  the  fact  that 
11©  C  II*.  Thus,  q(D)f  =  0. 

On  the  other  hand,  as  any  function  in  £(0),  /  can  be  written  in  the  form  /  = 
J2geeCde0'  hence  0  =  Q(D)f  =  Yloee  Q{6)c0eB-  Since  any  finite  set  of  exponentials  are 
linearly  independent,  we  conclude  that 

eg  ^0  =>  q(0)  =  0. 

Thus,  the  subset 

©'  :=  {6  e  ©  :  eg  ^  0} 

lies  in  the  hyperplane  q  =  0.  Ranging  the  directional  derivative  Dx  over  an  (m  —  2)- 
dimensional  space  (which  is  possible  since  p  is  bivariate)  we  obtain,  as  in  the  proof  of  (b), 
that  ©'  lies  on  a  circle.  It  remains  to  show  that  #©'  >  n  —  1.  However,  since  /  €  E(Q'), 
and  p  =  /j.  p  €  n©-,  and  therefore  its  space  of  derivatives  D(p)  lies  in  II©<.  as  well.  Since 
dimD(p)  was  shown  to  be  n  —  1,  we  conclude  that  #©'  =  dimn©'  >  dim  D(p)  =  n  —  1. 


5.  Strict  positive  definiteness  of  K  for  S1:  an  analytic  approach 

In  the  case  of  interpolating  on  the  circle,  dim'H^  =  2  for  all  A:  =  1,2,...,  and  H°k  is 
spanned  by  the  two  functions  cos  kr  and  sin  kr,  where  (r,  r)  are  the  polar  coordinates  in 
IR2.  These  well-known  facts  can  be  nicely  used  in  the  course  of  study  of  Problem  2.8.  We 
will  connect  our  problem  to  the  distribution  of  zeros  of  bounded  univariate  exponentials, 
and  use  the  obtained  characterization  to  derive  separate  necessary  and  sufficient  conditions 
for  the  induction  of  strict  positive  definiteness  by  K  for  Sl. 

Throughout  the  section,  we  will  make  an  essential  use  of  the  following  univariate 
exponential  space: 

£n-=  {  ^2  cTeiT  :  #T  <  n,  cT  €  IR,  Vr}. 

t€TC[0,2tt) 

Also,  £  :=  Un>  i £n.  Note  that,  importantly,  we  allow  only  real  coefficients  cT  in  the 
definition. 

We  first  record  and  prove  the  basic  observation  that  will  be  utilized  throughout  this 
section. 
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Theorem  5.1.  Given  K  C  7Z+,  the  following  conditions  are  equivalent. 

(a)  K  induces  strict  positive  definiteness  for  Sl  of  order  n. 

(b)  K  is  a  uniqueness  set  for  £n,  i.e.,  only  the  0  function  in  £n  vanishes  identically  on  K. 
Proof.  Selecting  the  basis 

y'/fc)(cos  r,  sin  r)  :=  cos  kr,  cos  r,  sinr)  :=  sin  kr 

for  Hi,  and  writing  a  given  finite  ©CS1  in  the  form  ©  =  {(cos  r,  sin  r)  :  r  €  T  C  [0.  27t) } , 
a  tj-pical  column  of  the  matrix  B k  of  (2.4)  then  has  the  form 

(cos  kr,  sin  kr)T . 

Therefore,  given  c  €  IRT,  we  conclude  that 

=  (Y  cr  coskr)'  +  (^2  CrS'mkr)2  =  |  ^  cTelkr |2. 

reT  r£T  r£T 

Introducing  the  exponential 

/:x4  ^22  cretTX, 
re  T 

we  conclude  that  \\Bkc\\^  =  |/(A)|2.  Now.  suppose  that  I\  induces  strict  positive  definite¬ 
ness  of  order  n,  and  /  is  any  exponential  of  above  form  in  £n.  With  (Bk)k  the  corresponding 
matrices,  we  have,  for  some  k  €  A',  Bkc  ^  0,  and  hence,  at  that  k  f(k)  /  0.  The  converse 
implication  is  obtained  by  reversing  the  argument.  4 

The  last  theorem  leads  to  the  following  necessary  condition  for  the  induction  of  strict 
positive  definiteness  by  I\  C  7L+. 

Corollary  5.2.  In  order  that  I\  C  7Z+  induce  strictly  positive  definite  for  Sl,  it  is  nec¬ 
essary  that  K  has  infinite  intersection  with  any  set  of  the  form  k7L+,  k  €  IN.  The  same 
applies  to  sets  of  the  form  k/2  -f  k7Z+,  provided  that  k  is  even. 

Proof.  In  order  to  unify  the  proof,  we  assume  that  the  relevant  arithmetic  progression 
is  of  the  form  Z  :=  a  +  VLL+.  a  <  k.  We  will  further  assume  that  K  has  only  finite 
intersection  with  Z ,  and  will  use  that  to  construct  a  linear  combination  f*  of  exponentials 
(eir)re[o,2rr)  that  vanishes  on  I\.  The  crux  in  the  proof  is  that,  if  a  =  0  or  a  =  k/2,  the 
coefficients  in  the  representation  /*  =  cTCir  are  real.  Thus,  for  these  cases,  /*  €  £, 
and  the  desired  result  is  then  implied  by  Theorem  5.1. 

Set 

T:=  {{2irl)/k  :  /  =  0, 1,  1}, 
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and  define  the  univariate  exponential 

/  :=  £eiT(-  -a). 
r6  t 

Then  f(l)  =  0,  for  0  <  l  <  k  -  1,  /  ^  a,  and  since  /  is  fc-periodic.  we  conclude  that  / 
actually  vanishes  on  7L+\Z. 

By  our  assumption,  the  set  K  H  Z  is  finite.  Let  n  denote  its  cardinality,  and  let 
i0,  t\.t2,  •  •  •  ,<2n  be  chosen  in  a  way  that  (a)  each  set  ti  +  T  is  a  subset  of  [0, 2ir).  and 

(b)  the  sets  ti  4-  T,  l  =  0, _ 2 n  are  pairwise  disjoint.  The  restriction  to  K  n  Z  of  the 

2 n  4-  1  exponentials  {eit(/}^0  must  be  linearly  dependent  over  the  reals,  since  the  set 
K  n  Z  contains  no  more  than  n  points.  Let  f*  be  a  non-trivial  linear  combination  with 
real  coefficients  of  these  functions  that  vanishes  on  K  n  Z .  Furthermore,  since  /  vanishes 
on  2Z+\Z.  so  does  every  eUlf .  and  therefore  /*  vanishes  on  that  set.  as  well.  This  implies 
that  f*  vanishes  on  K.  too.  Finally,  since  the  spectra  of  the  exponentials  {e^,/}/  are 
pairwise  disjoint,  f*  cannot  be  identically  0  (since  any  finitely  many  exponentials  are 
linearly  independent ) . 

It  remains  to  show  that  f*  G  £,  and  this  is  the  part  where  we  need  the  assumption 
(Y  =  (),A-/2.  Indeed,  for  such  choice  of  a,  we  observe  that,  the  numbers  e;T(-a),  r  G  T  are 
real,  hence,  in  these  cases,  in  /  G  £ ,  implying  that  e,t,/  6  £,  too.  Since  £  is  a  linear  space 
over  the  reals,  we  conclude  that  f*  lies  in  that  space,  too. 

Consequently,  we  have  found  a  noil-trivial  exponential  f*  €  £,  that  vanishes  on  A'.  4 

Note  that  the  exponential  /*  that  was  used  in  the  proof  of  the  above  result  lies  in 
£(2n+  i)k-  This  raises  the  following  question:  assuming  I\  to  contain  exactly  n  numbers 
from  the  set  a  +  k7Z+,  is  it  possible  that  K  induces  strict  positive  definiteness  on  Sl  of 
orders  smaller,  still  closer,  to  (2n  +  l)k  (or,  roughly  speaking,  is  the  construction  of  /*  in  the 
last  proof  uses  as  few  exponentials  as  possible)?  Some  answer  to  this  question  is  provided, 
for  the  set  1  +  22Z+,  by  the  choice  K  :=  {0, ...,2n}.  Here,  by  Theorem  3.2,  K  induces 
strict  positive  definiteness  of  order  An  +  1,  which,  in  view  of  Theorem  5.1,  implies  that  no 
non-zero  exponential  in  £^+1  can  vanish  on  K.  On  the  other  hand,  A  contains  n  numbers 
from  the  set  1  4-  22Z+.  Thus,  for  this  A',  we  have  used  (2n  +  1)2  =  An  4-  2  frequencies  to 
construct  /* ,  while  the  present  argument  shows  that  we  could  not  use  4n  4- 1  frequencies. 
In  summary,  at  least  for  that  example,  the  construction  of  /*  was  “economical” . 

Our  next  result  is  a  stronger,  more  quantitative  version,  of  statement  (c)  in  Theorem 
3.5. 

Theorem  5.3.  Let  K  be  a  subset  of  7L+  and  let  n  be  a  positive  integer.  Assume  that, 
for  every  .7  €  Cs\{l}  of  cardinality  n  —  1,  there  exist  a  G  ZZ+  and  k  G  IN\(  JIN)  such  that 
{a,  0  4-  Jfc, . . . ,  a  4-  («  -  1  )k}  C  A'.  Then  K  induces  strict  positive  definiteness  of  order  n 
on  Sl. 
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The  following  lemma  is  required  for  the  proof  of  the  theorem. 

Lemma  5.4.  Let  T  oe  any  finite  subset  o/(0, 2ir).  Then  the  set  of  integer  zeros  of  the 
function  F  :=  YlreT^^  ~  1)  JS  °f  the  form  JTL,  with  J  a  subset  of  IN\1  of  cardinality 
<#T. 

Proof.  Let  r  €  T.  The  condition  eiT(k)  =  1  is  equivalent  to  k  G  ~7L.  Since  ~  >  1, 
the  group  ZZ  fl  ^-TL  is  a  proper  subgroup  of  TL.  Thus,  the  zero  set  of  F  is  the  union  of  at 
most  #T  proper  subgroups  of  ZZ+.  4 

Proof  of  Theorem  5.3.  Assuming  K  to  satisfy  the  assumptions  made  in  the  theo¬ 
rem.  we  will  also  assume  that  there  exists  an  exponential  /  G  £n  that  vanishes  on  K.  and 
will  show  that  the  two  assumptions  contradict  each  other.  Our  claim  would  then  follow 
from  Theorem  5.1. 

Being  in  Sn,  the  exponential  /  has  the  form 

/=^creir,  T  C  [0, 27t),  ct  G  IR\0, 

rex 

with  #T  <  n.  Without  loss,  we  may  assume  that  0  €  T,  and  that  Co  =  1  (otherwise,  we 
divide  /  by  cTociTo,  with  r0  the  smallest  number  in  T).  By  Lemma  5.4,  the  positive  integer 
zeros  of  the  function 

F  II  («*  -  » 

r€T\0 

are  of  the  form  JIN,  with  J  C  1N\1  of  cardinality  <  n.  By  our  assumptions  on  K,  there 
exist  a  €  Z5+,  and  k  G  IN\(  JIN),  such  that  (a,  a  +  k, . . . ,  a  +  (n  —  l)fc}  C  K.  For  that  k, 
F{k)  ^  0,  that  is  1  £  {eir(A;)}reT\o. 

Let  p  be  a  polynomial 

p  :  1 1-4  ^  cljF 

j 

whose  zero  set  is  {e;r(fc)  :  r  G  T\0}.  Then,  degp  <  n,  and  p(l)  ^  0.  Let  p(V)  be  the 
induced  difference  operator 

p(V)  :  g  ^  o,jg(-  4-  jk). 

3 

Then,  p(V)ejT  =  ejr  p(ejr(fc))  =  0,  for  each  r  G  T\0.  Hence  p(V)/  =  p(l)  ^  0,  i.e.,  p(V)/ 
is  a  non-zero  constant.  However,  since  we  assume  that  /  vanishes  at  each  of  the  points 
ct  +  jk,  j  =  0, . . . ,  n  —  1,  it  follows  that  p(V)/  must  vanish  at  a.  Thus,  we  have  reached 
the  desired  contradiction.  4* 

The  value  k  =  1  is  always  admissible  in  Theorem  5.3.  Hence: 
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Corollary  5.5.  If  I\  C  Z£+  contains  n  consecutive  integers,  then  K  induces  strict  positive 
definiteness  of  order  n  on  S1. 


Another  special  case  is  the  following: 

Corollary  5.6.  Assume  that  A  =  K7Z+.  Then  it  induces  strict  positive  definiteness  on 
S1,  provided  that  it  is  not  of  the  form  J2Z+,  #J  <  oo. 

Proof.  We  invoke  Theorem  5.3.  For  that,  we  let  n  be  any  positive  integer,  and  J  be 
any  finite  subset  of  EV\{1 }.  By  our  assumption,  there  exist  k  G  K\{JJS).  The  property 
K  =  K7L+  therefore  implies  that  {k,  2k, . . . ,  nk}  c  K.  By  virtue  of  Theorem  5.3.  K 
induces  strict  positive  definiteness  of  order  n  for  S1.  Since  n  was  arbitrary,  our  claim  is 
proved.  4* 

Finally,  we  complement  Corollary  5.2  by  the  following  result  (which  completes  the 
proof  of  Theorem  3.5): 

Corollary  5.7.  The  set  I\  :=  2Z+\(n  +  k7L+)  induces  strict  positive  definiteness  for  Sl, 
provided  that  a  £  |2Z+. 

Proof.  Thanks  to  Theorem  5.1,  it  suffices  to  show  that  no  exponential  /  G  €  vanishes 
on  K. 

Let  /  =  ]Crg Tcreir  be  an  exponential  that  vanishes  on  I\.  We  will  show  that, 
necessarily,  one  of  the  coefficients  (cr)  is  not  a  real  number,  implying  thus  that  /  £  £ . 
Without  loss  (see  the  proof  of  Theorem  5.3)  we  can  assume  that  0  G  T,  and  that  co  =  1. 
Further,  we  can  also  assume  without  loss  that  To  :=  {0,2ir/fc, . . .  ,2n(k  -  l)/fc}  C  T 
(otherwise,  we  may  add  the  missing  exponentials  with  0  coefficients).  For  r  G  To,  and 
t  G  T,  the  difference  operator  Vt  :  g  ►  g  —  eit(-k)g(-  +  k )  annihilates  e u,  and  satisfies 
Vt(eiT)  =  (1  -  eit(—k))eiT  (since  ejr  is  ^-periodic).  Therefore,  applying  V  :=  IlfgTVTo^ 
to  /,  we  obtain,  with  c  :=  rii€T\To(l  ~  £«(— &))  7^  0, 


V/  =  c]T  cTeiT. 

t€  To 


On  the  other  hand,  it  is  clear  that  V/  still  vanishes  on  K,  and  in  particular  V/  vanishes 
on  (0,1,...,  k  —  1}\«.  But,  the  exponentials  {ejT  }t£T  are  the  characters  of  the  group 
{0, . . . ,  k  —  1}.  hence  they  are  linearly  independent  over  that  set,  and  consequently,  since 
c0  =  1,  cT  =  eir(-a),  r  G  To.  However,  our  assumption  on  a  and  k  ensures  that  ^ 
or,  equivalently,  that  c27r/fc  =  e2;Tt/fc(aO  is  not  a  real  number. 
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6.  Strict  positive  definiteness  of  K  for  Sm  J:  an  algebraic  approach 

Here,  we  attack  the  problem  from  a  completely  different  angle.  While  the  core  of  the 
argument  used  in  the  previous  section  was  the  connection  between  2-dimensional  harmonic 
polynomials  and  their  analytic  completions,  the  course  here  is  of  algebraic  nature.  It  is 
based  on  the  realization  of  harmonic  polynomials  as  the  kernel  of  the  laplacian,  thereby 
exploits  the  connection  between  a  homogeneous  polynomial  ideal  (here,  the  principal  one 
generated  by  T)  and  its  kernel  in  II  (here,  the  harmonic  polynomials).  We  will  require  here 
the  notations  which  were  introduced  before  Theorem  3.6.  Also,  recall  that  a  subspace  H 
of  A  or  II  is  homogeneous  if  each  of  the  maps  k,  k  £  ZZ+,  maps  H  into  itself. 

Some  of  the  results  in  this  section  can  be  developed  in  a  more  general  (and  in  our 
opinion  more  natural)  setup. 

Definition  6.1.  Let  Cl  be  a  subset  of  JRm,  H  a  homogeneous  subspace  of  II.  and  n  a 
positive  integer.  We  say  that  H  is  correct  of  order  n  on  Cl  if  for  any  ©  C  Cl  of 
cardinality  n,  and  any  function  F  defined  on  ©  there  exists  a  polynomial  p  £  H  that 
interpolates  F  (on  Q). 

To  see  the  connection  between  the  new  definition  and  our  original  problem,  one  chooses 
il  :=  Sm~l,  and  H  :=  Uk  (defined  as  in  the  introduction).  Then,  as  asserted  in  Theorem 
2.10,  the  induction  by  I\  of  strict  positive  definiteness  of  order  ti  on  5m~1  is  equivalent  to 
the  nth-order  correctness  of  H  on  5m_1. 

Theorem  6.2.  Let  H,  n  and  Cl  be  as  in  the  above  definition.  Then  H  is  correct  of  order 
n  on  12  if  and  only  if  for  every  exponential  f  £  En(Cl),  there  exists  k  £  7Z+  such  that  the 
kth  order  homogeneous  differential  operator  k(f)(D)  does  not  annihilate  H  D  11°. 

Proof.  For  9  £  IRm,  let  60  be  the  functional  60  :  f  •->  f(9),  and  for  a  finite  0  C  !Rm, 
let  <5©  be  span{^0}0€©.  H  fails  to  be  correct  of  order  n  on  u  if  and  only  if  for  some  0  C  Cl 
of  cardinality  n,  dim/!)©  <  dim <5©  =  n.  But  in  (and  only  in)  such  a  case,  there  would  be 
X  =  Ce^9  e  orthogonal  to  H ,  i.e., 

^cep(0)  =  0,  V  peH. 
ee© 

Defining  /  :=  c0e0,  we  obtain  an  exponential  /  €  En(Cl)  such  that  p(D)f( 0)  =  Xp  = 

0,  for  every  p  €  H.  For  p  £  H  D  n£,  we  have 

0  =  p(D)j(  0)  =  p(D)(k(f)m  =  (k(f)(D))p(  0)  =  (k(f)(D))p. 

Therefore,  the  condition  p(D)f( 0)  =  0,  all  p  €  f/,  is  equivalent  to  k(f)(D)  annihilating 
H  fl  11°  for  every  k.  ♦ 
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Proof  of  Theorem  3.6.  By  Theorem  2.10  (and  in  view  of  Definitions  2.8  and  6.1)  the 
required  property  of  K  is  equivalent  to  Hk  being  correct  of  order  n  for  Sm~l.  Therefore, 
in  view  of  Theorem  6.2.  we  need  to  prove  that  the  condition  stated  in  the  present  theorem 
is  equivalent  to  the  following:  "There  exists  /  €  En(Sm~l)  such  that,  for  every  k  €  K , 
( k(f))(D )  annihilates  'Hft.-nn®.  =  Since  H°k  is  the  kernel  in  11°  of  the  laplacian  T(D), 
the  last  condition  is  equivalent  to  the  divisibility  of  k(f)  by  T.  6 

In  the  rest  of  the  section,  we  prove  Theorem  3.7.  We  divide  the  proof  into  a  sequence 
of  several  lemmas,  some  of  them  might  appear  to  be  of  independent  interest.  Note  that 

CT-T2)  =dimn»_. 

Lemma  6.3.  The  operator  F(D)  induces  an  isomorphism  between  rn+1II  and  Fnn.  n  = 

0,1 . 

Proof.  First,  we  recall  that  II  is  the  direct  sum  of  the  spaces  {TkH)ke'Z+  ■  (A  quick 
proof  of  that  would  go  as  follows.  Since  H  is  the  kernel  in  II  of  the  homogeneous  differential 
operator  T(Z)),  II  is  the  direct  sum  of  H  and  the  principal  ideal  TII: 

n  =  h  e  rn. 

Multiplying  the  above  equation  by  Tk,  and  using  induction  we  obtain 

fc-l  k 

n  =  (Y,  ©  r*n  =  QT  rjn)  ©  r*+1n.) 

j=0  j= 0 

Next,  one  checks  directly  that,  for  each  k  €  ZL+,  there  exists  a  constant  cn,k  such  that 
the  operator  cn)A;rr(D)  is  the  identity  on  rn+1%°.  This  implies  that  r(D)  maps  Tn+1H 
1-1  onto  TnH. 

Finally,  the  decomposition  result  asserted  in  the  first  paragraph  of  the  proof  allows  us 
to  write 

rn+1n  =  ©*>  n+iTkH. 

The  desired  result  then  follows  from  an  application  of  T(D)  to  both  sides  of  the  last  equal¬ 
ity,  and  invoking  the  isomorphism  assertion  from  the  second  paragraph  of  the  proof. 

4 

Lemma  6.4.  Assume  that  f  €  A  and  satisfies  T(D)f  =  f.  Assume  further,  that  for  some 
k,  ti  £  2Z+ , 

(6.5)  (A-  +  2j)(f)  €  m,  j  =  0, . . . ,  n. 
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Then, 


(6.6)  (k  +  2n)(/)  €  rn+1n. 

Proof.  By  induction  on  n.  If  n  =  0,  the  claim  in  (6.6)  is  assumed  in  equation  (6.5). 
Assume  thus  that  the  claim  is  valid  for  n  —  1  >  0.  Since  T(D)f  =  /,  T(D)((k  4-  2 n)(/))  = 
(k  4-  2n  —  2)(/),  hence  by  the  induction  hypothesis 

r(D)((F+2^)(/))  e  rnn. 

Invoking  Lemma  6.3,  we  conclude  that 

(k  +  2n)(f)  eu  +  rn+1n. 

/ 

The  result  then  follows  from  the  assumption  that  ( k  +  2 n)(f)  is  divisible  by  T  (and  the 
fact  that  no  polynomial  in  'H\ 0  is  divisible  by  IT).  4 

Proof  of  Theorem  3.7.  Let  0  C  5m_1  of  cardinality  n  be  given,  and  let  /  E  E(Q). 
By  Theorem  3.6  (more  precisely,  by  the  statement  made  in  the  paragraph  that  follows 
Theorem  3.6),  we  need  to  show  that  K{f)  is  not  divisible  by  T.  Assume,  to  the  contrary, 
that  Tv  (/)  is  divisible  by  I\  YVe  first  note  that,  for  some  k  <  n/2,  say,  k0,  k0(f)  #  0. 
(Indeed,  by  the  definition  (4.2)  of  IIq,  if  k(f)  =  0  for  all  k  <  n/2,  then  ne  must  contain  a 
polynomial  of  degree  >  n/2,  in  contradiction  with  Lemma  3.1).  Since  /,  as  any  function  in 
E(Sm~l),  satisfies  T(D)f  —  f,  one  concludes  from  the  fact  that  k0(f )  ^  0,  that  i(f )  ^  0, 
for  every  i  e  k0  +  27L+.  Because  fe0  <  n/2,  our  assumptions  on  K  imply  that  there  exist 
j  integers  A:,  k  -  2,  k  -  4, . . . ,  k  -  2j  +  2  in  K,  that  further  satisfy  k  -  2 j  +  2  ~  ko  €  2ZZ+. 
Thus,  the  fact  that  i{f)  ^  0  holds,  in  particular,  for  any  integer  i  in  this  progression.  On 
the  other  hand,  each  of  these  integers  lies  in  K,  and  since  I\  (/)  6  I'll,  Lemma  6.4  implies 
that 

k(f)  e  Fn. 

As  the  rest  of  the  proof  will  establish,  this  last  conclusion  contradicts  the  fact  that  /  € 

E(Q). 

Since  /  is  a  linear  combination  of  {ee}ee©,  we  have  that 

£(/)  = 

oeT 

with  T  C  0,  and  with  any  9,6'  €  T  being  linearly  independent.  Here, 

(6-)k  :x^(6-  x)k. 
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Let  T'  be  a  non-spanning  set  in  X  of  maximal  cardinality.  Since  dimlT^  >  <r(0)  > 
#(X\X'),  there  exists  a  non-trivial  homogeneous  polynomial  pi  €  that  vanishes  on 
T\T'.  Writing  pi  in  the  form  pi  =  Tlp  with  p  £  HI  and  i  >  0,  we  conclude  that  p  is  a 
homogeneous  polynomial  of  degree  r  <  j  which  vanishes  on  T\T' .  An  application  of  p(D) 
to  k(f)  yields  a  linear  combination  of  the  form 

£  a„(0.)‘-r. 
eeT’ 

Since  V  does  not  span,  the  above  polynomial  is  of  <  m  variables. 

On  the  other  hand.  k(f)  =  rsP,  for  some  P  €  n\(m),  s  >  j.  Therefore,  for  some 
c  ^  0, 

P(D)(k(/))  =  P(D)(rp)  =  cPr‘~'p  +  r*-’'+1p,. 

Since  cpP  is  non-zero  and  is  not  divisible  by  T,  we  conclude  that  p(D)(k(f))  is  a  non-zero 
polynomial  in  m,  hence,  in  particular,  cannot  be  of  less  than  m  variables,  and  we  thus 
have  reached  the  desired  contradiction.  4 
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